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Abstract 



I. INTRODUCTION 
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We study numerically the influence of 
momentum-conserving dephasing on the trans- 
port in a disordered chain of scatterers. Loss of 
phase memory is caused by coupling the trans- 
port channels to dephasing reservoirs. In con- 
trast to previously used models, the dephas- 
ing reservoirs are linked to the transport chan- 
nels between the scatterers, and momentum con- 
serving dephasing can be investigated. Our 
setup provides a model for nanosystems exhibit- 
ing conductance quantization at higher temper- 
atures in spite of the presence of phononic inter- 
action. We are able to confirm numerically some 
theoretical predictions. 
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The recently discovered jumps in the con- 
ductance of metallic nanowires and nanocon- 
tacts |I| have demonstrated impressively 
that quantum effects can dominate transport 
properties of small structures even at ambi- 
ent temperature. This is generally explained 
by the lateral confinement which induces a 
rather large subband spacing of the trans- 
port modes of the conductor 0. Within a 
single transport mode, phononic backscatter- 
ing is considerably reduced because of the 
large amount of momentum transfer of J9 = 
2hkp that has to be provided by a phonon. 
However, there remains the question whether 
the underlying picture of "ideal conductors" 
0,^ is appropriate to describe e.g. metallic 
nanocontacts; if one insisted in a description 
in terms of ideal waveguides, one should ac- 
count for the unavoidable imperfections of 
real systems like impurities or geometric devi- 
ations from the ideal-waveguides setup. The 
inclusion of coherent scattering centers in cor- 
responding transport models leads, however, 
to the almost immediate loss of quantization 
of the conductance in the absence of some 
stabilizing mechanism. 

An interesting question is whether de- 
phasing, induced by phononic interaction, 
could provide such a mechanism: Interac- 
tion with longitudinal phonons, although not 
causing momentum transfer, can well be at 
the origin of dephasing which counteracts the 
weak localization of coherent backscattering. 
In this way, temperature could be helpful to 



re-establish quantization in this specific sit- 
uation. For an investigation of this question 
transport models are needed which comprise 
elastic scattering as well as dephasing, while 
allowing for vanishing momentum transfer. 

Many approaches to dephasing are based 
on reservoirs which destroy phase informa- 
tion in the same way as momentum informa- 
tion of passing electrons [^j-Q ; therefore, cor- 
responding transport models are not suited 
for the investigation of the above discussed 
transport regime. 

For the present study we have developed 
a model where coherent and incoherent scat- 
tering processes are conceptually separated. 
Coherent scattering is described by general 
elastic scatterers. We investigate single-mode 
transport through a one-dimensional chain of 
scatterers (Fig. |l|). 

Dephasing is induced using a common ap- 
proach employing virtual electron reservoirs 
^. We have linked the reservoirs to the 
transport channels between elastic scatterers, 
exploiting the fact that electrons in trans- 
port channels automatically possess definite 
momenta. This seems to be the easiest way 
to implement momentum-selective coupling 
to a reservoir as has also been discussed by 
Datta [10|. Chemical potentials of such reser- 



voirs thus correspond to occupation numbers 
of a local momentum distribution function as 
it occurs in the Boltzmann equation. This 
implies furthermore that emission from the 
reservoirs is uncorrelated with respect to the 
opposite direction of motion. 

As depicted in Fig. |I|, electrons absorbed 
by one of the two reservoirs between two scat- 
terers are fed back incoherently into either 
the same transport channel or into the ad- 
jacent channel. In the latter case, they con- 
serve their momentum as in the case of co- 
herent propagation, while having, however, 
lost their phase memory. This permits us 
to investigate continuously different degrees 
of momentum conservation, from the case of 



full momentum conservation to the case of 
momentum flip, independently from the de- 
phasing strength itself. 

Since only two scattering channels per 
scatterer are involved, our model is also ef- 
ficient for numerical calculations which is 
of practical interest in order to cope with 
ensemble-averaging. It can also be easily ex- 
tended towards multi-mode transport. 

In the following section, we briefly gen- 
eralize a method for the calculation of the 



scattering matrix of a composed system |Tl 



towards the inclusion of dephasing reservoirs 
linked to "inner" channels. Section III re- 
views the involved characteristic lengths for 
phase, momentum and coherent localization 
and their relationship to the parameters of 
our model. Numerical results for the metal- 
lic and insulating regime are presented in 
section IV, and the influence of momentum- 
conserving dephasing on the transport is dis- 
cussed. 



II. NUMERICAL APPROACH 

We consider a one-dimensional chain con- 
sisting of N scatterers (see Fig. p. Each scat- 
terer has two channels and is described by 
its scattering matrix S^ which relates linearly 
the outgoing to the incoming amplitudes at 
the scatterer. We have used the parametriza- 
tion 
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where r^ denotes the reflection probability 
of the scatterer k, k = 1, . . . ,N. A sys- 
tem of non-interacting scatterers would be 
described by a 2N x 2N scattering matrix 
K containing the S^ on its diagonal. De- 
noting the 2 A^-tuples of incoming and out- 
going amplitudes by c+ and c_, respectively, 
we have c_ = Kc+. For linked scatterers, 
we may distinguish between external and in- 
ternal channels; external channels are con- 



nected to the outside, e.g., the left and the 
right contact in the present case (channels 1 
and 2N), while internal channels are inter- 
connected: Amplitudes outgoing into inter- 
nal channels are propagated to incoming am- 
plitudes in internal channels, thus acquiring 
a phase shift exp(ip) which may comprise ei- 
ther a simple geometric phase p = qd for a 
wave vector q and a distance d, or, generally 
also a magnetic phase j A ■ ds in presence of 
a vector potential A. This can be described 
by an operator P which propagates internally 
outgoing amplitudes while annihilating exter- 
nally outgoing amplitudes, c+_int = Pc_. Us- 
ing the decomposition c+ = c+_int + c+,ext, the 
relation 
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is readily obtained ^^. The wanted scatter- 
ing matrix S of the composed system relates 
incoming to outgoing amplitudes in external 
channels, i.e., 
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it is thus given by 
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where the indices i, j label only the external 
channels. 

The scattering matrix S then yields the 
transmission matrix T according to Tij = 
ISijl"^; Tij is the transmission probability be- 
tween the external channels j and i. The 
transmission matrix thus relates tuples of 
(dimensionless) outgoing and incoming cur- 
rents. 
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The associated incoming and outgoing cur- 
rents are given by the the squared moduli of 
the amplitudes in the external channels. 

To incorporate dephasing in this ap- 
proach, we assign a virtual reservoir to each 



internal channel / and denote by 7^ the prob- 
ability of getting absorbed by such a reser- 
voir. The current into this reservoir is thus 
given by 7; |c_ p. Current conservation then 
demands a corresponding loss term being ac- 
counted for also by the operator P; therefore, 
the assigned phase factor for coherent prop- 
agation of c_ has now to be multiplied by 
a/1 — 7;. For the considered chain it is given 
by e.g. 



Pmi = exp{ipmi)\/l -li 
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where rra = / ± 1 for / even or / odd, respec- 
tively, pmi denoting the phase for coherent 
propagation. Speaking in terms of currents, 
only a fraction 1 — 7; of the current outgoing 
in channel / is propagated coherently while a 
fraction 7^ is absorbed by the reservoir. In the 
following, we assume 7^ = 7 throughout for 
all internal channels / while we can formally 
assign 7„ = 1 to external channels n which 
are fully absorbing; in this way, the virtual 
reservoirs are treated in the same way as the 
external contacts. 

In the presence of the virtual reservoirs, 
the dimension of the transmission matrix 
T has become larger and the matrix ele- 
ments describing transmission from the vir- 
tual reservoirs are also needed. They can 
be obtained by replacing c+^ext in Eq. ^ with 
c+,int^/l where c+^jnt is given by the l-th. unit 
vector for virtual reservoir /. The squared 
moduli of c_ (multiplied by 7 for internal 
channels i) then yield the elements Tn, i.e., 
column / of T. 

The virtual reservoirs can now be used 
to introduce dephasing into the system [Q. 
For this purpose, electrons absorbed by the 
virtual reservoirs have to be re-emitted inco- 
herently. The motion of electrons inside the 
considered system (see Fig. |I]) can then be de- 
scribed as a Markovian process (see e.g. |0]): 
They are injected by the left contact, pass 
through virtual reservoirs and finally get ab- 
sorbed by either the right or the left contact. 



Adopting this point of view, we may renounce 
the notation of a chemical potential f^ which 
is convenient since we are interested in the 
total transmission only. 

It is at the level of the transmission ma- 
trix where current conservation is imposed. 
We describe momentum conservation by a 
parameter a ranging from to 1; electrons 
absorbed by a reservoir between two scat- 
terers in the chain conserve their momen- 
tum with probability a while their momen- 
tum gets flipped with probability a = 1 — a 
(see Fig. |ID. The total transmission T is then 
given by 
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as sum of transmission probabilities from the 
left to the right contact with 0, 1, ... inter- 
mediate passages through virtual reservoirs. 
This infinite series describing the random 
walk of the electrons can be summed analyti- 
cally. We note that T itself is symmetric due 
to time reversal invariance. Therefore, it is 
no surprise that the virtual reservoirs flip the 
momenta of absorbed electrons if the latter 
are fed back into the same channel; incoming 
and outgoing amplitudes in the same channel 
correspond to opposite momenta. Thus, our 
procedure of imposing current conservation 
only onto a pair of channels is well beyond 
the widely used description of Biittiker 0; 
it is, however, covered by the more general 
approach of Datta 



III. CHARACTERISTIC LENGTHS 

The phase coherence length lq> is given by 
the mean distance that an electron travels 



coherently, i.e., the distance it conserves its 
phase. In order to relate it to the absorption 
probability 7, we consider the unperturbed 
system (without disorder, i.e., vanishing elas- 
tic backscattering) and obtain /$ in units of 
the number of passed dephasing regions be- 
tween the scatterers as 



U 



i:n=i<i-ir 



E 



71=1^ 



7)" 



1 

7' 



The mean free path l^ is the mean distance 
that an electron travels without being scat- 
tered, i.e., without momentum change. In or- 
der to obtain Im in the unperturbed system, 
we consider the (incoherent) reflection prob- 
ability f of a single dephasing region which 
is given by f = 07. The resulting (dimen- 
sionless) resistance caused by a pair of corre- 
sponding reservoirs is po = ''^/(l ~ ''^)i which 
is additive for two incoherenty coupled scat- 
terers 1^. We now show that Im '■= Po^ is a 
useful definition for the mean free path: For 
a chain of A^ dephasing areas, this yields a 
serial (four-probe) resistance of p = Nl 
and, because of p = (1 — T)/T, a total trans- 
mission of 
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Equation ^ describes an Ohmic length depen- 
dence of the transmission; a section of length 
Im contributes a resistance unit to the total 
resistance 



10 



The effect of the coherent disorder can be 
described by the localization length l^. We 
are extracting this length from the ensemble- 
averaged transmission of purely coherent sys- 
tems, using the quantum scaling law of An- 
derson et al. (Ref. fl^). 
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which differs from the well-known Landauer 
result by a factor of 2. Ref. [|T^ predicts, 
assuming random phases between the scat- 
terers, a localization length of the ensemble. 



measured in number of scatterers, as 
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with the ensemble averages of the reflection 
and transmission probabilities r and t = 1 — r 
of a single scatterer. For our numerical inves- 
tigation, we have employed a constant reflec- 
tion probability Vk = r and random phases 
P2k+i,2k = P2k,2k+i e [0, 27r[ (indicated by arcs 
in Fig. 0). 

Performing numerical tests by logarith- 
mically averaging up to 10^ realizations of a 
purely coherent system (7 = 0) consisting of 
up to A^ = 5 X 10"^ scatterers, we have a found 
quite a good agreement with Eqs. |T^, |TT| (see 
Fig. H). We have also tried a random box dis- 
tribution with average r instead of constant 
reflection probability; this neither had an in- 
fluence on the validity of Eq. 0, nor on the 
following results involving both coherent and 
incoherent processes. 



IV. RESULTS 

In order to investigate the interplay be- 
tween the "microscopic" parameters r, 7 and 
a and their counterparts, the characteristic 
lengths /^, /$ and Im, we distinguish between 
the "conducting" regime /$ < /^ and the "in- 
sulating" regime where /^ < /$. 

Two representative results are shown in 
Fig. ^ (left and right, respectively), where the 
computed ensemble-averaged total transmis- 
sion is plotted as a function of the number 
of scatterers for dephasing which conserves, 
randomizes and flips momenta {a = 1, 0.5 
and 0, respectively). Transmission of the 
fully coherent system (7 = 0) is similar to 
Fig. 1^ (but not exactly the same due to the 
smaller number of averaged realizations) . 

The "conductor" can be thought of be- 
ing composed of phase-coherent units which 
are linked incoherently, resulting in an Ohmic 



length dependence of transmission. Since de- 
phasing destroys localization, dephasing with 
full momentum conservation (a = 1) always 
enhances transmission when compared to a 
fully coherent system (7 = in Fig. ^. 

The widely discussed case of simultane- 
ous phase and momentum randomization 
corresponds to a = 0.5; here, with increas- 
ing number of scatterers, transmission first 
drops due to incoherent backscattering, but 
there is always a certain length A^ such that 
the enhancement of transmission due to sup- 
pression of coherent localization outperforms 
the inelastic backscattering for systems larger 
than A^. The absolute value of the trans- 
mission, however, already has become very 
small. Therefore, enhanced conductance due 
to dephasing can only be expected for a large 
degree of momentum conservation a. 

In the insulating regime, the dephasing 
rate, but also the backscattering rate are 
much smaller than in above discussed, "con- 
ducting" case. With l^ < /$, conduction 
properties are dominated by localization ef- 
fects and cannot be understood in terms of in- 
coherently coupled phase coherent units. De- 
phasing supports transmission significantly 
only for system sizes larger than /$ when a 
deviation from the exponential dependence 
(described by Eq. [Tl|) in the right diagram 
of Fig. 1^ can be noticed. Finally, one would 
also expect an Ohmic behaviour for very large 
system sizes A^ ^ /$. Some enhancement 
of transmission in the insulating regime com- 
pared to the conducting regime appears but 
only for a < 1 (which can be seen comparing 
the two plots in Fig. 3). This is not unphysi- 
cal: A large /$ means small dephasing, which 
leads to a suppression of the quantum local- 
ization of the coherent system and thus in- 
creases transmission. For even larger dephas- 
ing (approaching the conducting regime), we 
get an increasing additional Ohmic resistance 
which finally causes the transmission to drop 
again (for a < 1). 



In the insulator, the relative enhancement 
of the transmission for a = 1 compared to 
a = 0.5 or a = is much smaller than for 
the case of the conductor, momentum con- 
servation has less influence in this regime. 

The importance of momentum conserva- 
tion is best seen in Fig. ^ where the depen- 
dence of the transmission on the dephasing 
rate 7 is plotted. Only for a = 1, transmis- 
sion increases monotonously with increasing 
7. In the incoherent limit 7 = 1, an Ohmic 
transport regime is reached with each single 
scatterer contributing the resistance l7^ (see 
Eqs. ^, |TT]), resulting in a transmission 

T = 1/(1 + N/l^ + N/l^). (12) 

For a = 0.5 and a = 0, transmission passes 
through a certain maximum which has al- 
ready been found for the case a = 0.5 in 
Ref. 0; beyond this maximum, the addi- 
tional resistance of the dephasing reservoirs 
causes a rapid decline. 

In the following, we compare our numeri- 
cal results (Figs. ^ and ^ to a formula which 
has been proposed by Band et al. |l^ on the 
basis of ensemble averaged quantities. The 
dashed lines in Figs. |^ and ^ have been cal- 
culated using a slightly modified version of 
this formula. 

In this approach, the transmission is given 
as 
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with pd being the resistance due to disor- 
der for the case of infinite /„• Equation [13 
can be seen as a generalization of Eq. ^ to- 
wards the inclusion of disorder. The validity 
of Mathiessen's law is assumed in the sense 
that the two processes which cause momen- 
tum change, namely momentum fiip by the 
reservoirs and coherent refiection by the scat- 
terers, give independent contributions to the 
resistance. Note that they do not depend on 
l^ and a, respectively. Surely, quantum resis- 
tance is qenerally not additive. However, the 



range of validity of the approach of Ref. 13 is 
surprisingly large (see Fig. 4). Our numerical 
verification of the result of Band et al. shows 
excellent agreement. A pre- requisite for this 
is our dephasing model, which allows for mo- 
mentum conservation. 

In the presence of dephasing (for finite 
^$^)) Pd can be built up recursively from 
the resistance of the purely coherent case, 
Pc{x;l^) = e^/'« — 1, and from the probability 
density P{x; N,l^,l^) for an electron to travel 
coherently a distance x within the disordered 
wire and to lose phase information at posi- 
tion X by a dephasing event. Both pc{x]l^) 
and P(x; A^, /g, /#) are ensemble averages. In 
the case /$ < A^, Band et al. obtain 
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Since both localization and dephasing pro- 
cesses lead to an exponential length decay of 
the probability. 
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with X G [0,A^], is a plausible choice, with 
a disorder modified phase coherence length 
l\{l^,l<^). In their original work. Band 
et al. have proposed 
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however, in more recent work |T^,|T5[ , a phase 
coherence length in a disordered quantum 
wire is defined in a similar context as 
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Our numerical results are better described 
using Eq. |l^. 

Equation |I3| has been employed for Figs. |^ 
and ^. This gives a rather good agreement, 
except for the case l^ > l^ (right diagram 
of Fig. 1^) where the contribution of inelas- 
tic backscattering {a = 0) to the resistance 



is slighly underestimated; Mathiessen's law 
may thus not be valid in this regime. 

For the interesting case of enhanced de- 
phasing rates, Eq. |l^ can be used to estimate 
the achievable change of transmission com- 
pared to purely coherent transmission. The 
maximum enhancement AT is obtained for 
/$ = OJm = cxd; using Eqs. |l^ and |l^, we 
obtain with x '■= ^/h 

AT = l/(l + x)-e-^. (18) 

This yields a maximum ATmax ~ 0.20 for 

V. CONCLUSIONS 

We have presented a model of general 
scatterers which allows us to investigate 
quantum transport through a disordered sys- 
tem in the presence of dephasing with a tun- 
able degree of momentum conservation. We 
obtained the rather obvious result that mo- 
mentum conservation is essential if dephasing 
is to enhance transmission. From Fig. ^, it is 
evident that support for a quantized conduc- 
tance by dephasing can be expected in the 
regime Iq, < N k.1^. 

Our numerical results have been com- 
pared to a general but heuristic picture 
developed in Ref. [Q which is based on 
the validity of Mathiessen's rule and an ad 
hoc assumption for a probability density 
P{x] N, l^, lq>)] it is not derived from a micro- 
scopic theory, but gained its plausibility from 
its validity in various limiting cases. Further- 
more, it is not evident that average quantities 
as Pc{x) and P{x; N, l^, lq>) can be employed 
in a recursive formula for the transmission, 
since sample averaging is not a simple com- 
mutative operation and since there are differ- 
ent types of averaging procedures which of- 



However, while the previously discussed 
approach by Band et al. yields a good agree- 
ment for the ensemble-averaged transmission, 
it cannot predict fluctuations. For the inves- 
tigation of conductance quantization effects, 
it is highly desirable to describe a single sys- 
tem or at least the applicability of results for 
ensemble averages to a particular realization. 
Here, our general quasi-microscopic model 
provides a helpful tool. In nanosystems 
exhibiting conductance quantization, elastic 
backscattering is reduced (but not absent) 
because of the larger subband spacing which 
is itself induced through the geometrical con- 
finement. At higher temperatures, one has 
to account for phononic interaction, notably 
with logitudinal phonons, leading to dephas- 
ing but not inelastic backscattering. As men- 
tioned in the introduction, one can describe 
this situation by momentum-conserving de- 
phasing which was not included in previous 
models. While we found an excellent agree- 
ment with the approach of Ref. 13 over a wide 
range, the domain of applicability to (quasi- 
) quantization in nanosystems would be for 
A^ < /# < /^. Work along these lines is 
in progress. An advantage compared to an 
heuristic approach is that we are able to in- 
vestigate also fluctuations and single realiza- 
tions. 



ten proved intuition wrong in the past [12 



Therefore, it is rather surprising that Eq. [T^ 
gives a good flt of our data, except for the 
strongly insulating regime, /^ ^ /$ ^ A^. 
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FIG. 1. A chain consisting of A^ = 3 scatter- 
ers 5i , ^2 and ^3 . Channels 1 and 6 link the sys- 
tem to the left and right contact. Electrons out- 
going from the scatterers in inner channels 2-5 
may travel coherently or incoherently between 
the scatterers. When traveling coherently, they 
aquire a phase shift p (indicated by arcs) . Alter- 
natively, they may get absorbed (dotted arrows, 
probability 7) between the scatterers by reser- 
voirs (wavy lines). From the reservoir they are 
re-emitted (dashed arrows) either with same mo- 
mentum (probability q) or inversed momentum 
(probability a = 1 — a). 
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FIG. 2. Transmission of a chain of N co- 
herently coupled scatterers with backscattering 
probability r = 0.01. Each point is obtained by 
logarithmically averaging 10^ realizations. The 
dashed line corresponds to Eq. ^ with a local- 
ization length /^ = 99 = (1 — r)/r. 
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FIG. 3. Transmission as function of the 
number of scatterers for three different de- 
grees of momentum conservation a. The elas- 
tic backscattering probabiUty is r = 0.01 corre- 
sponding to a locahzation length l^ = 99. The 
left diagram is calculated for dephasing proba- 
bility 7 = 0.05, the right diagram for 7 = 0.005. 
Each point is obtained by logarithmically aver- 
aging 1000 realizations. Also indicated is the 
transmission of the coherent system (7 = 0). 
The dashed and dotted lines indicate solutions 



of Eq. 13 and Eq. 11, respectively 
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FIG. 4. Total transmission for a chain of 
N = 500 scatterers as a function of the dephas- 
ing parameter 7 for three different degrees of mo- 
mentum conservation a. The elastic backscat- 
tering probability is r = 0.01 corresponding to 
a localization length l^ = 99; each point is ob- 
tained by logarithmically averaging 1000 realiza- 
tions. The lines are obtained from Eq. |l^ which 
is applicable for 7 > 1/A^ (indicated by the ar- 
row). 



